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ABSTRACT 

Through  the  Optimum  Path  Aircraft  Routing  System  (OPARS) , 
Fleet  Numerical  Oceanography  Center  (FNOC)  has  committed  it- 
self to  providing  a  computerized  flight  planning  service 
remotely  accessible  via  dial-up  communications  lines.   The 
question  arises  as  to  whether  the  proposed  number  of  tele- 
phone lines  will  be  adequate  to  provide  a  level  of  service 
previously  provided  by  the  Lockheed  Jetplan  system.   This 
study  provides  a  detailed  analysis  of  the  response  delays  for 
the  OPARS  flight  plan  system.   In  addition,  estimates  are 
given  of  communication  requirements  when  various  levels  of 
demand  prevail,  and  under  conditions  in  which  the  FNOC  computer 
is  busy  or  idle. 
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I.   SUMMARY 

In  the  proposed  OPARS  flight  plan  system,  users  at 
37  remote  terminals  (Figure  1)  connect  directly  to  the 
Fleet  Numerical  Oceanography  Center  (FNOC)  computer  system 
via  11  telephone  lines.   The  total  amount  of  time  that  a  line 
is  busy,  the  total  OPARS  service  time,  can  be  modeled  as 
the  sum  of  three  sub-system  service  times: 

1.  An  interactive  program  setup  time  T(setup); 

2.  A  delay  that  the  OPARS  program  experiences  in 
the  FNOC  computer's  input  queue  T( input  queue); 

3.  An  OPARS  program  run  time  T(run). 

The  expected  values  of  these  individual  sub-system  service 
times  can  be  computed  and  then  summed  to  provide  a  mean 
or  expected  service  time  for  the  entire  OPARS  request  as 
follows : 

E[T(service) ]  =  E[T(setup) ]  +  E[T(input  queue)]  +  E[T(run)] 
From  this,  Erlang's  formula  equation  (1)  can  be  employed  to 
compute   individual  long-run  state  probabilities  (the 
probability  that  i  lines  are  busy)^  p.,  and  ultimately  the 
probability,  p,.. ,  of  having  all  eleven  lines  busy. 

In  support  of  the  analytic  model,  a  computer  program 
was  written  to  simulate  the  entire  OPARS  flight  plan  request 
process,  from  initial  dial-up  through  program  completion. 
The  simulation  also  includes  the  arrival  and  servicing  of  other 
FNOC  computer  programs  which  interfere  with  the  processing 
of  the  OPARS  program. 
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Figure  1.   Remote  Terminal  Sites 
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Both  the  analytic  model  and  the  simulation,  computed 
state  probabilities  for  OPARS  demand  rates  of  2,  4,  6,  ...20 
per  hour  and  under  conditions  in  which  the  FNOC  computer  is 
busy  or  idle.   The  results  clearly  indicate  that  the  11 
telephone  lines  can  handle  demand  rates  up  to  three  times 
as  great  as  those  currently  being  experienced  by  the  Lockheed 
Jetplan  system  before  p   ,  the  probability  of  all  lines 
busy,  exceeds,  .05. 
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II.   INTRODUCTION 

Fleet  Numerical  Oceanography  Center  (FNOC)  is  currently 
testing  and  evaluating  a  computerized  flight  plan  system, 
referred  to,  for  short,  as  OPARS.   This  sytem,  developed  to 
replace  the  Lockheed  Jetplan  flight  plan  sytem,  provides  users 
at  remote  sites  with  direct  access  to  the  FNOC  computer 
via  11  telephone  lines.   The  purpose  of  this  study  is  to 
determine  if  the  intended  number  of  telephone  lines  would 
be  adequate  to  ensure  a  low  probability  of  having  all  lines 
busy. 

The  number  of  lines  busy  at  any  time  t,  {N.( t )  |  t>o}  can  be 

modeled  as  a  Birth-Death  process  with  inter-arrival  times 

that  are  assumed  to  be  independent  and  exponentially 

distributed  but  with  service  times  that  clearly  are  not. 

Fortunately,  in  a  communications  system  characterized  by 

calls  that  arrive  in  a  stationary  or  time-homogenous  Poisson 

Process  of  rate  A  and  vie  for  n  lines  with  no  queue  (blocked 

calls  lost),  the  long  run  probability  of  having  i  lines  busy, 

P-j_,  can  be  computed  from  Erlang's  formula. 

i 
n   -  1  (At)  /il  0<i<N  (1) 

pi  ~\ N  .  _. 


I  Ut)/j! 

j=0 


i>N 


This  formula  has  the  surprising  feature  of  being  valid  for 
any  service  time  distribution.   Indeed,  given  the  above 
assumption,  one  need  only  obtain  the  mean  service  time,  t, 
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in  order  to  calculate  p. .   The  majority  of  the  effort  in 
this  study  was  to  characterize  the  OPARS  system  in  such  a 
way  so  that  a  mean  service  time  could  be  derived  under  various 
operating  conditions. 

This  paper  begins  with  some  introductory  material 
concerning  the  FNOC  computer  systems, the  OPARS  flight  plan 
system  and  historical  usage  of  the  Lockheed  Jetplan  system. 
Section  IV  is  a  detailed  description  of  the  analytic  model 
used  and  a  discussion  of  the  supporting  siumulation  program. 
Finally,  Section  V  contains  the  tabled  results  and  an  analysis 

Notation  used  in  this  paper  is  consistent  with  that  found 
in  Queueing  Theory  literature  and  will  be  explained  whenever 
introduced.   In  addition,  a  listing  and  discription  of 
all  notation  and  abbreviations  can  be  found  on  page  seven. 
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III.   BACKGROUND 


A.   FNOC  COMPUTER  SYSTEM 

FNOC  currently  has  four  computers  in  service  at  their 
Monterey  facility,  but  only  one,  referred  to  as  HAL  is 
accessible  to  the  remote  OP ARS  terminals.   HAL  is  a  CDC-6600 
computer  with  two  central  processors  and  330  K  of  octal 
central  memory.   HAL  is  accessed  within  FNOC  by  approximately 
twenty-five  interactive  terminals  as  well  as  by  a  batch 
job  system  that  routinely  processes  hundreds  of  programs  a 
day.   Since  HAL  alone  is  involved  in  providing  service  to 
OPARS  remote  terminals,  we  will  only  consider  it  in  the 
discussion. 

1.   HAL  Batch  System 

HAL's  batch  system  is  composed  of  a  priority  ranked 
"first-in,  first-out"  input  queue,  which  can  be  considered 
exterior  to  the  computer,  and  a  priority  ranked  execute 
queue  from  which  jobs  are  selected  by  the  scheduler  for 
processing  (Figure  2).   When  a  job  enters  the  system,  it 
first  enters  the  input  queue,  at  which  the  job's  user-assigned 
priority  establishes  its  initial  position.   While  in  the 
queue  the  job  slowly  accrues  additional  "wait  time" 
priority  which  ensures  that  even  the  lowest  priority  jobs 
are  eventually  run.  The  largest  number  of  daily  jobs  are 
restricted  to  priorities  of  1,  2  or  3  (3  being  the  highest) 
but  priorities  up  to  77  can  be  assigned.   The  priority  of 
OPARS  jobs  is  fixed  automatically  at  60  which  immediately  puts 
it  ahead  of  all  but  a  few  other  jobs  in  the  queue. 
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When  space  is  available  within  the  computer  the  scheduler 
removes  the  first  job  (with  respect  to  priority  then  to 
time-of-arrival)  from  the  input  queue  and  moves  it  to  the 
execute  queue.   Upon  entering  the  execute  queue  the  job 
loses  all  of  its  previously  gained "wait-time"  priority  but 
immediately  begins  to  gain  it  back  as  the  job  waits  to  be 
processed.   In  this  queue  the  jobs  with  the  higher  assigned 
priorities  gain  this  additional  priority  faster  than  those 
with  low  assigned  priority.   This  mechanism  results  in  higher 
priority  jobs  getting  selected  for  processing  sooner  than 
low  priority  jobs.   When  selected  by  the  scheduling  routine, 
the  program  is  moved  into  central  memory  and  it  is  processed 
for  a  unit  length  of  time  or  until  it  attempts  to  access 
a  device  (disk,  tape,  etc.)  which  is  unavailable.   When  this 
event  occurs  the  job  is  removed  form  central  memroy  and 
returned  to  the  execute  queue,  having  its  accured  priority 
reset  to  zero.   The  process  continues  in  this  manner  until 
the  program  is  completed  and  it  is  transferred  to  the  out- 
put queue. 

2.   Current  FNOC  Workload 

From  available  central  processor  and  central  memory 
utilization  profiles  (Figure  3)  it  can  be  readily  seen  that 
resource  demands  on  HAL  can  be  separated  into  two  states : 
A  Busy  or  High  Demand  state  in  which  the  computer  is  operating 
at  maximum  capacity,  and  an  Idle  or  Low  Demand  state  in  which 
most  of  the  resources  are  immediately  available.   The  High 
Demand  state  is  characterized  by  a  full  execute  queue,  and  a 
lengthy  input  queue,  while  the  Low  Demand  state  exhibits  an 
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execute  queue  with  only  an  occasional  job  in  it  and  an  in- 
put queue  which  is  empty.   Fortunately  the  transition 
period  between  these  states  is  quite  brief  and  therefore 
a  transition  state  is  not  needed. 

The  High  Demand  period  coincides.,  not  surprisingly > 
with  normal  working  hours  except  that  it  extends  well  into 
the  evening,  often  as  late  as  midnight.   Any  jobs  entering 
the  system  during  this  period  must  be  delayed  for  a  time 
in  the  input  queue  before  being  run.   During  the  idle  period, 
on  the  other  hand,  jobs  pause  only  momentarily  in  the  input 
queue  before  being  run. 

B.   OPARS  PLIGHT  PLAN  SYSTEM 

The  OPARS  flight  plan  system  provides  the  user  with 
an  optimum  (with  respect  to  time)  route  of  flight  given  forecast 

weather  and  user  inputs  such  as  aircraft  type  and,  take  off 
weight.   For  flights  within  areas  that  have  a  defined 
routing  structure  the  optimization  is  fairly  easy  because 
only  a  few  routes  are  candidates  for  the  optimum.   Over 
water  or  point  to  point  flight  plans  are  more  difficult 
because  there  is,  in  theory,  an  uncountably  infinite  number 
of  available  routes. 

There  are  two  major  sections  of  the  OPARS  program:   an 
interactive  portion  and  the  optimization  program.   When  a 
user  at  a  remote  site  requests  an  OPARS  flight  plan  one 
first  sets  up  the  program  with  the  Interactive  Response 
Generator  (IRG).   Using  a  query  and  response  technique,  the 
IRG  prepares  the  main  program  by  obtaining  required 
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information  from  the  user.   During  this  session  the  IRG 
is  not  checking  input  information  for  validity,   but  only 
for  format.   For  example,  an  entry  of  ABCE,  as  the  four- 
letter  identification  code  for  the  destination  airfield, 
would  be  accepted  even  though  no  such  airfield  code 
exists.   When  all  of  the  required  information  is  entered, 
the  IRG  allows  the  user  to  review  and  change  any  input 
if  necessary.   Then,  upon  command ^  the  program  is  put  into 
the  HAL  computer  batch  system. 

As  discussed  in  the  previous  section,  the  OPARS  job 
enters  the  input  queue  along  with  all  other  FNOC  jobs  and 
awaits  its  turn.  Hereafter  the  OPARS  program  is  handled 
as  any  other  batch  program  with  one  exception.   Upon 
completion  the  flight  plan  is  returned  automatically  to  the 
terminal  where  the  request  was  entered.   If  for  some 
reason  the  line  has  been  disconnected^  the  flight  plan  is 
placed  into  an  output  file  from  which  it  can  be  retrieved 
later. 

C.   HISTORICAL  USAGE/PROJECTED  DEMAND 

In  an  effort  to  evaluate  expected  demand  for  the  OPARS 
flight  plan,  several  methods  were  employed.   They  were: 

1)  A  linear  regression  model  of  historical  usage  by 
the  Navy  to  estimate  demand  through  1981; 

2)  expectations  of  remote  site  users; 

3)  Lockheed's  records  of  previous  use. 

None  of  these  methods  alone  provided  the  required  level  of 
accuracy  and,  only  by  combining  the  three,  could  a  reasonable 
estimate  be  made. 
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The  Navy's  own  records  provided  the  best  overall 
information.   This  information,  shown  in  Figure  4,  consists 
of  monthly  totals  of  Navy  requests  for  the  Jetplan.   The 
trend  is  clearly  an  increase  in  usage  and  a  simple  linear 
regression  yielded  the  projected  demands  shown  in  Figure  5. 
These  monthly  totals,  unfortunately  did  not  provide  any 
information  as  to  how  requests  varied  during  the  day. 


1977 
Month(per  day) 

1978 

Month(per  day) 

1979 
Month(per  day) 

JAN 

866(28.6) 

1177(38) 

2033(65.6) 

FEB 

1012(3^.9) 

1107(39-5) 

1926(68.8) 

MAR 

1076(34.7) 

1856(59.9) 

2107(68) 

APR 

1169(33.9) 

1319(44) 

1887X62. 9) 

MAY 

1103(35.6) 

1490(48.1) 

2119(68.4) 

JUN 

1174(39.1) 

1224(40.8) 

2081(69.4) 

JUL 

1711(55.2) 

1146(37) 

2011(64.9) 

AUG 

1066(34.4) 

1365(44) 

2247(72.5) 

SEP 

977(32.3) 

1226(40.9) 

2218(73.9) 

OCT 

1142(36.8) 

1415(45.6) 

NOV 

1103(36.8) 

1321(44) 

DEC 

1107(35-7) 

1386(44.7) 

Figure  4.  Historical  Usage  of  Lockheed's  Jetplan 


Y  INTERCEPT 

28.79 

JUL  80  ESTIMATE  80 . 4  /DAY 

SLOPE  1.2 

JUL  81  ESTIMATE  94.8  /DAY 

COR.  COEFF. 

.826 

Figure  5.  Linear  Regression  of  Usage  Data 
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In  order  to  obtain  the  distribution  of  daily  requests, 
expected  demand  information  was  requested  from  all  remote 
site  users.   The  resulting  information  suffered  from 
imprecision  and  a  lack  of  completeness  (only  about  half  of 
the  remote  sites  responded) . 

As  a  final  effort,  an  attempt  was  made  to  secure  actual 
daily  records  from  the  Lockheed  Corporation.   The  Lockheed 
personnel  were  unable  to  provide  me  with  actual  data,  but 
over  the  course  of  several  interviews  a  picture  of  the 
daily  demand  distribution  was  pieced  together. 

The  results  of  these  data  collection  efforts  appear  in 
Figure  7.   These  graphs  show  the  projected  demand  for 
July  1980  and  July  1981  distributed  using  the  composite 
results  of  Lockheed  information  and  user  expectations. 
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IV.   THE  MODEL 


We  start  the  modeling  discussion  by  reviewing 
some  basic  assumptions  about  the  system. 

First  and  foremost  is  the  assumption  concerning 
time-homogeneous  Poisson  arrivals  for  OPARS  requests. 
Without  this  assumption,  the  Erlang  formula  (1)  would  not 
be  valid  and  a  much  more  complex  model  would  be  required. 

Secondly,  the  service  times  for  the  three  sub-systems 
are  assumed  to  be  mutually  independent. 

Finally,  as  discussed  in  the  last  section,  the 
conditions  for  the  computers  workload  will  be  divided  in- 
to two  distinct  periods.   These  are  the  High  Demand  period 
with  a  nonempty  input  queue,  and  the  Low  Demand  period 
where  the  OPARS  requests  bypass  the  input  queue  and  go 
directly  into  the  computer. 

A.   THE  ANALYTIC  MODEL 

Since  the  main  thrust  of  the  analytic  effort  is  to 
obtain  an  overall  mean  service  time,  an  immediate 
simplification  would  be  to  somehow  break  the  OPARS  system 
into  two  or  more  sub-systems  whose  service  times  would 
be  more  easily  obtained.   With  this  in  mind,  we  start  by 
separating  the  service  into  the  interactive  sub-system 
and  the  batch  sub-system.   The  latter  is  still  rather 
formidable  because  of  the  nature  of  the  input  and 
execute  queues  and  the  complexity  of  the  scheduling 
routine.   The  final  step  is  to  reduce  the  batch  system 


into  an  input  queue  sub-system  and  an  execution  sub-system. 
The  reasons  for  this  may  not  be  readily  apparent  as  there 
are  other  modeling  techniques  such  as  an  M|G|m  multiserver 
queueing  system  with  priorities,  which  would  handle  the 
system  as  a  whole.   Such  a  system  would  require  that  each 
job  priority  have  a  mean  service  time.   Unfortunately,  there 
is  no  correlation  between  a  jobs  priority  and  its  core 
and  central  processor  requirements.   Additional  unrealistic 
assumptions  would  be  required  which  would,  in  the  end, 
reduce  the  validity  of  the  result. 

1.  The  Interactive  Sub-System 

First  to  be  considered  is  the  interactive  sub-system. 
An  estimate  for  the  mean  time  for  this  sub-system  is 
readily  comput  ble  from  the  available,  albeit  scarce,  data. 
This  data(Figure  7)  was  obtained  by  timing  FNOC  technicians 
on  trial  OPARS  runs.   The  mean  of  these  data  is  5-1  minutes 
and  will  be  assumed  to  be  the  same  during  both  High  and 
Low  Demand  periods. 

2 .  Input  Queue  Sub-System 

The  next  sub-system  to  model  is  the  input  queue  sub- 
system.  Under  low  demand  conditions  the  mean  service  time 
is,  by  definition,  equal  to  zero.   During  the  High  Demand 
period  it  is  not  quite  so  simple.   Again}  by  previous 
definition,  the  High  Demand  period  has  a  continuously  non- 
empty input  queue.   This  very  simple  assumption  suggests  that 
the  input  queue  may  be  modeled  as  a  single-server  queueing 
system  by  itself,  the  server  being  the  first  position  in  the 
queue.   One  must  merely  postulate  the  interarrival  process 
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41911 

7'22" 

6r2" 
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5'09" 

5T17" 

5T01" 

4»2i" 

3?57" 

415511      . 

X 

=5-1      MIN 

a 

=    1.0 

Figure  7.  Interactive  Setup  Times 

and  the  service  process,  the  latter  representing  the  times 
between  successive  departures  from  the  input  queue.   If  we 
can  now  model  the  arrival  rate  as  Poisson  and  the  service 
distribution  as  exponential  we  have  an  M/M/l  queueing 
system  ranked  by  priority  for  which  there  is   a  simple 
closed  form  solution  for  expected  waiting  times.   The 
interested  reader  is  urged  to  consult  Griffin   [2]  for  the 
derivation  of  this  result. 

There  was  no  data  available  to  determine  the  precise 
manner  in  which  OPARS  jobs  enter  the  input  queue  and  so  a 
Poisson  arrival  rate  is  as  plausible  as  any  other.   However 
the  service  distribution  (the  interdeparture  times  from  the 
input  queue)  can  be  measured  by  observing  a  real  time  display 
of  HAL's  input  queue  and  timing  departures.   The  resulting 
data  (Figure  8)  and  histogram  (Figure  9)  are  surprisingly 
close  to  the  exponential  distribution.   In  fact,  the  Chi- 
squared  goodness  of  fit  test  resulted  in  a  test  statistic  of 
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Figure  8.  Input  Queue  Interdeparture  Time 
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Figure  10.   Higher  Priority  Job  Counts 


27 


0.533 5  which  indicates  that  the  data  is  very  close  to  being 
exponential  with  AQ  =  53-52  per  hour. 

The  final  consideration  is  the  priority   ranking. 
As  stated  earlier  an  OPARS  job  is  assigned  a  priority  of  60 
while  FNOC  jobs  can  have  priorities  from  1  to  77.   However, 
the  batch  jobs  of  lower  priority  can  be   ignored  and  we  are 
left  with  only  higher  priority  jobs  to  consider.   Available 
data  (Figure  10)  showed  higher  priority  jobs  arriving  at  the 
rate  of  3-^5  per  hour.   The  one  remaining  assumption  is  that 
these  higher  priority  jobs  also  arrive  according  to  a 
Poisson  Process. 

The  expression  for  OPARS  expected  input  queue  is 
then, 
E[T( INPUT  QUEUE)]  =  ■ ~ ■ _  (2) 

y 


y  (1  -  *n,+  *H  (  1  -  XhJ 


where  u  is  the  service  rate,  u  is  the  OPARS  arrival  rate  and 
X        is  the  higher  priority  job  arrival   rate.   This  expression 

a 

is  valid  only  when: 

XH  +  Xy    >  1. 

V 
This  input  queue  waiting  time  was  computed  for  OPARS  arrival 

rates  of  2,  4,  6,..., 20  per  hour  with  results  listed  in 

column  three  of  Figure  16. 


3.   The  Execution  Sub-Systems 

In  the  third  and  final  section,  the  execution 
sub-system,  we  will  characterize  the  OPARS  program  run  time  as 
well  as  delays  due  to  other  factors. 
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Because  of  the  central  memory  requirements  of  the  OPARS 
program  (150K)  only  two  OPARS  jobs  are  allowed  in  the  computer 
at  one  time.   Any  additional  requests  arriving  during  such 
a  state  would  be  required  to  wait  in  the  input  queue  and  would 
be  passed  over  by  the  scheduling  routine  as  it  selected  jobs 
for  processing.   Because  the  two  OPARS  jobs  can  be  processed 
simultaneously,  a  queueing  system  of  the  form  GI/G/2  is 
suggested.   This  type  of  system  is  uncomfortable  to  work 
with,  so  simplifying  assumptions  will  be  made. 

First,  because  of  a  lack  of  data  to  the  contrary, 
we  will  again  ssume  Poisson  arrivals  of  OPARS  jobs.   It 
now  remains  to  describe  the  service  distribution.   Unfortunately, 
even  with  Poisson  arrivals,  multi-server  systems  with  any  but 
exponential   service  times,  are  difficult  to  analyze  math- 
ematically. 

Figures  11  and  12  contain  actual  OPARS  run-time  data 
for  Low  Demand  and  High  Demand  periods  respectively.   These 
data  reflect  the  time  an  OPARS  job  spends  in  the  system 
from  it's  transfer  to  the  execute  queue  until  it's  completion. 
Histograms   of  these  data  are  in  Figures  13  and  14.   and  are 
clearly  not  exponential  .   However,  if  the  data  can  be 
considered  to  be  from  an  Erlang  distribution  then  the  nomo- 
graph in  Figure  15  (see  Hillier  and  Leiberman  [4])  can  be 
used  to  obtain  an  approximate  value  for  N,  the  mean  number 
of  OPARS  jobs  in  the  system.   From  this  the  mean  service  time 
can  be  computed  via  Little's  result: 

E[T(run)]  =  5  (3) 

o 
where  ^0is,  as  before,  the  arrival  rate  for  OPARS  jobs. 
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Figure  11  .  OPARS  Program  Run  Times:  Low  Demand 
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Figure  12.   OP'ARS  Program  Run  Times:  High  Demand 
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Figure  15.  Nomograph  of  the  Long  Run  Average 
Number  of  Customers,  N,  in  an  M/E  /2  Queue 
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Figure  16  .    Busy   System  Summary    (Time    in   Minutes) 


Xo 

E[T( setup)] 

E[T 

2/HR 

5.1 

0 

4/HR 

5.1 

0 

6/HR 

5.1 

0 

8/HR 

5.1 

0 

10/HR 

5.1 

0 

12/HR 

5.1 

0 

14/HR 

5.1 

0 

16/HR 

5-1 

0 

18/HR 

5.1 

0 

20/HR 

5.1 

0 

E[T(input  queue)]     p  /  N  /  E[T(run)]     E[T( Services)] 


074/. 15/4. 5 

9.6 

148/. 30/4.5 

9.6 

222/. 46/4.6 

9.7 

296/. 63/4. 7 

9.8 

369/. 83/5.0 

10.1 

443/1.1/5.5 

10/6 

517/1-4/6.0 

11.1 

591/1.7/6.4 

11.5 

665/2.1/7.0 

12.1 

738/2.8/8.4 

13.5 

Figure     17.    Idle   System  Summary    (Time    in  Minutes) 


The  run  data  was  parameterized  as  an  Erlan^ 
distribution  with  the  following  results : 


BUSY 

IDLE 

r 

4 

3 

A 

.6645 

.6772 

1   r 

6.02  min 

4.43  min 

Figure  18.   Erlang  Distribution  Parameters 
Where  —  is  the  mean  run  time. 

^R 

Finally,  the  expected  execution  service  time, 
E[T(runs)],  was  computed  for  a  series  of  OPARS  arrivals, 
as  before,  with  the  results  in  column  4  of  Figures  16-  and  17. 
4.   Summary 

In  summary  we  have  now  succeeded  in  representing 
the  total  OPARS  service  time  as  the  sum  of  three,  readily 
computable  sub-system  times  (Figure  19) : 

E[T(service) ]  =  E[T(setup)]  +  E[T(input  queue)]  +E[T(run)1.  (4) 
Erlang* s  formula  can  now  be  employed  to  compute  expected 
state  probabilities  for  lines  busy. 
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B.   THE  SIMULATION 

As  an  alternative  solution  technique,  two  computer 
programs  were  written  to  simulate  the  OPARS  system,  one  each 
for  High  and  Low  Demand  computer  states.   Output  of  the 
simulations  were,  as  in  the  analytic  model,  state  probabilities 
for  lines  busy  under  various  demand  states. 

There  are  two  key  differences  in  the  logic  flow  for  the 
two  simulations.   First  OPARS  requests  arriving  during  the 
busy  state  are  sent  to  the  input  queue  where  they  must  wait 
to  be  served.   Under  idle  conditions  jobs  go  directly  to 
the  computer.   Secondly,  in  the  busy  state,  higher  priority 
jobs  are  being  created  which  interfere  with  the  processing 
of  OPARS  jobs.   The  following  is  a  list  of  the  various 
distributions  and  their  parametric  values  used  in  the 
simulation. 

OPARS  Interarrival  time:  EXP  (x=  2 , 4, 6 , . . . 20/HR) 

High  Priority  Jobs  Interarrival  time:..  EXP(x=  3.45/HR) 

IRG   service  times: NORMAL  (y=  5.1  Min. ,  a=  1) 

Input  queue  Interdeparture  times :EXP(A=  53-5/HR) 

OPARS  run  times,  Idle:    GAMMA  U=  .721,  r  =  4.3*0 

OPARS  run  times.,  Busy:    GAMMA  (X=  .701,  r  =  3-D 

Figure  20:   Probability  Distribution  Used  in  Simulation 
Program 

The  simulations  were  run  for  48  hours  for  each  OPARS 
demand  rate  and  under  each  computer  state. 
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The  only  significant  factor  included  in  the  simulations 
but  not  in  the  analytic  model  had  to  do  with  a  retrial 
population.   In  the  event  of  all  lines  busy,  the  analytic 
model  assumes  arriving  requests  disappear, but  in  reality, 
people,  when  faced  with  a  busy  signal,  hang  up  and  try 
again.   The  simulation  retries  all  balked  attempts  after 
ten  minutes. 
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V.   RESULTS  AND  ANALYSIS 

Erlang  formula  (1)  state  probabilities,  pi,  were 
computed  for  OPARS  request  rates  of  2, 4, 6,..., 20  per  hour 
under  High  and  Low  Demand  conditions.   Appendix  A  contains 
the  FORTRAN  program  used  to  calculate  these  state  probabilities 
given  demand  rate  and  mean  service  time,  E[T( service) ] . 
Similarly  Appendices  D  and  C  contain  the  SIMSCRIPT  programs 
which  simulate  the  OPARS  system. 

The  results  of  both  methods  are  in  Figures  2.  and  22. 
The  state  probabilities  are  suprisingly  close  (differing 
by  less  than  2%   in  most  cases)  and  clearly  show  that  with 
eleven  operating  telephone  lines,  demand  rates  must  be  far 
in  excess  of  projected  requirements  (Section  III.C),  before 
the  probability  of  all  lines  being  busy  is  of  concern. 

One  should  be  aware,  however,  that  in  spite  of  the 
closeness  of  the  results  these  two  solutions  merely  serve 
to  verify  each  other.  Because  the  OPARS  system  is  not  fully 
functional  at  the  time  of  this  report,  there  is   no  way 
to  validate  either  method  against  the  actual  system. 

In  addition,  if  any  significant  changes  are  made  to 
either  the  OPARS  system  or  to  the  FNOC  computer  system, 
a  new  evaluation  would  be  required. 
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APPENDIX  A 
COMPUTER   PROGRAM   TO    CALCULATE   ERLANG   FORMULA   STATE   PROBABILITIES 


DIMENSION  PROB  (12,10),  Serve  (10) 

ZL  =   .0001 

NN  =  1 
5   READ  (5,10)  (SERVE (I), I  =  1,10) 
10   FORMAT (IP10. 5) 

LINE  =  12 

DO  60  J  =  1,10 

RHO  =  2.0  *  J/SERVE(J) 

DENOM  =0.0 

DO  30  K  =  1,LINE 

FACT  =1.0 

DO  20  L  =  1,K 

FACT  =  FACT  *  (L-l) 

IF  (L  .EQ.l)  FACT  =1.0 
20   CONTINUE 

DO  50  I  =  1,LINE 

FACTNU  =1.0 

DO  40  L  =  1.1 

FACTNU  =  FACTNU  *  (L-l) 

IF(L  .EG.  1)  FACTNU  =  1.0 
40   CONTINUE 

PROB  (I, J)  =  ( RHO** (1-1) /FACTNU) /DENOM 

IF  (PROB(I,J)  .LT.  ZL)  PROB(I,J)  =  0.0 
50   CONTINUE 
60   CONTINUE 

WRITE  (6,90) 

DO  70  I  =  1,  LINE 

M  =  1-1 

WRITE(6,65)  M,  (PROB(I,J),  J  =  1,10) 
65   FORMAT  (2X,I2,10F11.6) 
70   CONTINUE 

IF  (NN.EQ.2)  STOP 

NN  =  2 

GO  TO  5 
90   FORMAT  (Tl) 

END 
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APPENDIX  B 

SIMULATION  PROGRAM:  HIGH  DEMAND  STATE 

PREAMBLE 

NORMALLY  MODE  IS  INTEGER 

EVENT  NOTICES  INCLUDE  HIGH. PRI .JOB ,  OPARS.JOB,  RETRY, 

INPUT.  QUEUE, 

TEMPORARY  ENTITIES 

EVERY  JOB  HAS  A  TIME. OP. ENTRY  AND  A  PRIORITY  AND  MAY 

BELONG  TO  THE  QUEUE 

DEFINE  QUEUE  AS  A  FIFO  SET  RANKED  BY  PRIORITY 

THE  SYSTEM  OWNS  THE  QUEUE 

DEFINE  IRG,  THRUPUT,  QUEUE. TIME  AND  TIME. OF  ENTRY  AS 

REAL  VARIABLES 

DEFINE  NO. REQUESTS,  NO. ATTEMPTS,  PRIORITY,  LINES. 

BUSY,  RUN,  NO. BATCH,  SYSTEM,  BUSY  AND  EMPTY  AS  VARIABLES 

DEFINE  SUMMARY  AS  A  2-DIMENSIONAL,  REAL  ARRAY 

ACCUMULATE  STATE. PROB(0  to  11  BY  1)  AS  THE  HISTOGRAM  OF 

LINES, BUSY 

ACCUMULATE  QUEUE .STATS (0  TO  11  BY  1)  AS  THE  HISTOGRAM 

OF  N. QUEUE 

TALLY  AVE. QUEUE. TIME   AS  THE  AVERAGE  OF  QUEUE. TIME 


MAIN 


RESERVE  SUMMARY(*,*)  AS  12  BY  10 

LET  RUN  =  1 

LET  BUSY  =  1 

SCHEDULE  AN  OPARS.JOB  NOW 

SCHEDULE  A  STATISTICS  IN   2880  MINUTES 

START  SIMULATION 


END 


EVENT  OPARS.JOB 

SCHEDULE  AN  OPARS.JOB  IN  EXPONENTIAL. F( 30 . /RUN, 4)  MINUTES 

ADD  1  TO  NO. ATTEMPTS 

IF  LINES. BUSY   EQUALS  11 

SCHEDULE  A  RETRY  IN  10  MINUTES 

RETURN 

OTHERWISE 

ADD  1  TO  LINES. BUSY 

ADD  1  TO  NO. REQUESTS 

LET  IRG  =  NORMAL.F(5.1,1-0,5) 

SCHEDULE  AN  INPUT. QUEUE  IN  IRG  MINUTES 

RETURN 
END 

EVENT  HIGH. PRI. JCB 

SCHEDULE  A  HIGH. PRI. JOB  IN  EXPONENTIAL. F( 17. 4, 3)  MINUTES 

CREATE  A  JOB 

LET  PRIORITY  (JOB)  =  2 

FILE  JOB  IN  QUEUE 

RETURN 
END 
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EVENT  RETRY 

IF  LINES.  BUSY  EQUALS  11 

ADD  1  to  NO . ATTEMPTS 

SCHEDULE  A  RETRY  IN  10  MINUTES 

RETURN 

OTHERWISE 

ADD  1  TO  NO. REQUEST 

ADD  1  TO  LINES. BUSY 

LET  IRG  =  NORMA. F(5. 1,1. 0,5) 

SCHEDULE  AN  INPUT. QUEUE  INIRG  MINUTES 

RETURN 
END 

EVENT  INPUT. QUEUE 

IF  SYSTEM  EQUALS  EMPTY  AND  NO. BATCH  IS  LESS  THAN  2 

SCHEDULE  A  JOB . COMPLETION  IN  GAMMA. F( 4 . 43, 3 . 1 ,7)MINUTES 

ADD  1  TO  NO. BATCH 

OTHERWISE 

CREATE  A  JOB 

LET  TIME. OF. ENTRY (JOB)  =  TIME.V 

LET  PRIORITY(JOB)  =  1 

FILE  JOB  IN  QUEUE 

REGARDLESS 

RETURN 
END 

EVENT  EXECUTION 

IF  SYSTEM  EQUALS  BUSY 

SCHEDULE  AN  EXECUTUION  IN  EXPONENTIAL. F ( 1. 12,2)  MINUTES 

REGARDLESS 

IF  N. QUEUE  EQUALS  0 

RETURN 

OTHERWISE 

IF  PRIORITY(F. QUEUE)  =  2 

REMOVE  FIRST  JOB  FROM  QUEUE 

DESTROY  THE  JOB 

RETURN 

OTHERWISE 

IF  NO .BATCH  .  EQUALS  2 

RETURN 

OTHERWISE 

REMOVE  FIRST  JOB  FROM  QUEUE 

LET  QUEUE. TIME=  (TIME.V  -TIME. OF. ENTRY (JOB)  *  1440. 

DESTROY  THE  JOB 

IF  SYSTEM  EQUALS  BUSY 

SCHEDULE  A  JOB. COMPLETION  IN  GAMMA. F( 6. 06 , 4 . 34 , 6) MINUTES 

OTHERWISE 

SCHEDULE  A  JOB . COMPLETION  IN  GAMMA. F( 4 . 43 , 3 . 1, 7) MINUTES 

REGARDLESS 

ADD  1  TO  NO. BATCH 

RETURN 
END 
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EVENT  JOB. COMPLETION 

SUBTRACT  1  FROM  LINES. BUSY 

SUBTRACT  1  FROM  NO. BATCH 

IF  N. QUEUE  EQUALS  0 

RETURN 

OTHERWISE 

IF  SYSTEM  EQUALS  BUSY 

CANCEL  THE  EXECUTION 

REGARDLESS 

SCHEDULE  AN  EXECUTION  NOW 

RETURN 
END 

EVENT  STATISTICS 

START  NEW  PAGE 

PRINT  1  LINE  WITH  RUN*2  THUS 
ARRIVAL  RATE  =  **HOUR 

SKIP  3  OUTPUT  LINES 

PRINT  1  LINE  THUS 

PR(BUSY  LINES)    PR( INPUT  QUEUE) 

SKIP  2  OUTPUT  LINES 

FOR  I  =  1  TO  12,  DO 

PRINT  1  LINE  WITH  1-1,  STATE . PROB( I)/2  AND  QUEUE. 

STATS ( I) /2  THUS 
**       *****        ***** 

LOOP 

SKIP  5  OUTPUT  LINES 

PRINT  3  LINES  WITH  AVE. QUEUE. TIME, NO. ATTEMPTS  AND  NO. 

REQUESTS  THUS 
AVERAGE  OPARS  QUEUE  TIME  =  **.**** 
NO.OPARS   ATTEMPTED    **** 
NO. OPARS   COMPLETE  **** 

FOR  I  =  1  to  12,  DO 

LET  SUMMARY (I, RUN)  =  STATE. PR03( I) /2 

LOOP 

SCHEDULE  A  STATISTICS  IN  2880  MINUTES 

RESET  TOTALS  OF  LINES. BUSY,  N. QUEUE  AND  QUEUE. TIME 

IF  RUN  EQUALS  10 

GO  TO  FINAL. STATS 

OTHERWISE 

ADD  1  TO  RUN 

LET  NO. ATTEMPTS  '  =  0 

LET  NO. REQUESTS   =  0 

RETURN 

'FINAL. STATS' 

START  NEW  PAGE 

PRINT    5    LINES   THUS 
EMPTY   SYSTEM   SUMMARY 
STATE     2/HR  4/HR  6/HR     8/HR     10/HR     12/HR  14/HR     16/HR     18/HR     20/HR 

FOR  I  =  1  to  12,  DO 

PRINT  1  LINE  WITH  1-1,  SUMMARY  (1,1),  SUMMARY  (1,2), 

SUMMARY (1,3),  SUMMARY  (1,4), SUMMARY (1,5), SUMMARY (1,6), 

SUMMARY (I, 7),  SUMMARY  ( I , 8) , SUMMARY ( I ,9) , SUMMARY ( I ,10) ,THUS 
**   ****   ****   ****   ****   ****   ****   ****   ****   ****  •  **** 

LOOP 
STOP 
END  46 


APPENDIX  C 

SIMULATION  PROGRAM:  LOW  DEMAND 

PREAMBLE 

NORMALLY  MODE  IS  INTEGER 

EVENT  NOTICES  INCLUDE  HIGH. PRT. JOB ,OPARS .JOB, RETRY, INPUT. 

QUEUE, 

TEMPORARY  ENTITIES 

EVERY  JOB  HAS  A  TIME. OF. ENTRY  AND  A  PRIORITY  AND  MAY 

BELONG  TO  THE  QUEUE 

DEFINE  QUEUE  AS  A  FIFO  SET  RANKED  BY  PRIORITY 

THE  SYSTEM  OWNS  THE  QUEUE 

DEFINE  IRG,THRUPUT, QUEUE. TIME  AND  TIME. OF. ENTRY  AS 

REAL  VARIABLES 

DEFINE  NO  .  REQUESTS  ,'  NO.  ATTEMPTS  ,  PRIORITY,  LINES.  BUSY, 

RUN,  NO. BATCH,  SYSTEM,  BUSY  AND  EMPTY  AS  VARIABLES 

DEFINE  SUMMARY  AS  A  2-DIMENSIONAL,  REAL  ARRAY 

ACCUMULATE  STATE. PROB(0  TO  11  BY  1)  AS  THE  HISTOGRAM 

OF  LINES. BUSY 

ACCUMULATE  QUEUE. STATS .( 0  TO  11  BY  1)  AS  THE  HISTOGRAM 

OF  N. QUEUE  TALLY  AVE. QUEUE .TIME  AS  THE  AVERAGE  OF  QUEUE. 

TIME 
END 

MAIN 

RESERVE  SUMMARYC*,*)  AS  12  BY  10 

LET  RUN  =  1 

LET  SYSTEM  =  1 

LET  BUSY  =  1 

SCHEDULE  AN  OPARS.JOB  NOW 

SCHEDULE  AN  EXECUTION  IN  EXPONENTIAL. F ( 1, 12, 2) MINUTES 

SCHEDULE  A  HIGH. PRI. JOB  IN  EXPONENTIAL. F( 17 • 4 , 3)MINUTES 

CREATE  A  JOB 

LET  PRIORITY  (JOB)  =  2 

FILE  JOB  IN  QUEUE 

SCHEDULE  A  STATISTICS  IN  2880  MINUTES 

START  SIMULATION 
END 

EVENT  OPARS.JOB 

SCHEDULE  AND  OPARS.JOB  IN  EXPONENTIAL. F( 30. /RUN, 4)  MINUTES 

ADD  1  TO  NO. ATTEMPTS 

IF  LINES. BUSY  EQUALS  11 

SCHEDULE  A  RETRY  IN  10  MINUTES 

RETURN 

OTHERWISE 

ADD  1  TO  LINES. BUSY 

ADD  1  TO  NO.  REQUESTS 

LET  IRG  =  NORMA. F(5. 1,1. 0,5) 

SCHEDULE  AN  INPUT. QUEUE  IN  IRG  MINUTES 

RETURN 
END 

EVENT  HIGH. PRI. JOB 

SCHEDULE  A  HIG.PRI.JOB  IN  EXPONENTIAL. F( 17 . 4 , 3) MINUTES 
CREATE  A  JOB 
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APPENDIX  C 
SIMULATION  PROGRAM:  LOW  DEMAND  STATE  CONT'D 

LET  PRIORITY (JOB)  =  2 
FILE  JOB  IN  QUEUE 
RETURN 
END 

EVENT  RETRY 

IF  LINES. BUSY  EQUALS  11 

ADD  1  TO  NO  .ATTEMPTS 

SCHEDULE  A  RETRY  IN  10  MINUTES 

RETURN 

OTHERWISE 

ADD  1  TO  NO. REQUESTS 

ADD  1  TO  LINES  .BUSY 

LET  IRG  =  NORMAL. F(5. 1,1. 0,5) 

SCHEDULE  AN  INPUT. QUEUE  IN  IRG  MINUTES 

RETURN 
END 

EVENT  INPUT. QUEUE 

IF  SYSTEM  EQUALS  EMPTY  AND  NO. BATCH  IS  LESS  THAN  2 

SCHEDULE  A  JOB. COMPLETION  IN  GAMMA. F( 4 . 43 , 3 . 1 , 7)  MINUTES 

ADD  1  TO  NO  .BATCH 

OTHERWISE 

CREATE  A  JOB 

LET  TIME. OF. ENTRY (JOB)  =  TIME.V 

LET  PRIORITY(JOB)  =  1 

FILE  JOB  IN  QUEUE 

REGARDLESS 

RETURN 
END 

EVENT  EXECUTION 

IF  SYSTEM  EQUALS  BUSY 

SCHEDULE  AN  EXECUTION  IN  EXPONENTIAL. F( 1 . 12 , 2)  MINUTES 

REGARDLESS 

IF  N. QUEUE  EQUALS  C 

RETURN 

OTHERWISE 

IF  PRIORITY  (F. QUEUE)  -  2 

REMOVE  FIRST  JOB  FROM  QUEUE 

DESTORY  THE  JOB 

RETURN 

OTHERWISE 

IF  NO. BATCH  EQUALS  2 

RETURN 

OTHERWISE 

REMOVE  FIRST  JOB  FROM  QUEUE 

LET  QUEUE. TIME  =  (TIME.V  -  TIME. OF. ENTRY ( JOB) )  *  1440. 

DESTORY  THE  JOB 

IF  SYSTEM  EQUALS  BUSY 

SCHEDULE  A  JOB . COMPLETION  IN  GAMMA. F( 6 . 06, 4. 34, 6) MINUTES 
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SIMULATION  PROGRAM:  LOW  DEMAND  STATE  CONT'D 

OTHERWISE 

SCHEDULE  A  JOB . COMPLETION  IN  GAMMA. P( 4. 43,3 .1,7)  MINUTES 
REGARDLESS 
ADD  1  TO  NO. BATCH  ■ 
RETURN 
END 

EVENT  JOB. COMPLETION 

SUBTRACT  1FROM  LINES. BUSY 

SUBTRACT  1  FROM  NO. BATCH 

IF  N. QUEUE  EQUALS  0 

RETURN 

OTHERWISE 

IF  SYSTEM  EQUALS  BUSY 

CANCEL  THE  EXECUTION 

REGARDLESS 

SCHEDULE  AND  EXECUTION  NOW 

RETURN 
END 

EVENT  STATISTICS 

START  NEW  PAGE 

PRINT  1  LINE  WITH  RUN* 2  THUS 
ARRIVAL  RATE  =  **/HOUR 

SKIP  3  OUTPUT  LINES 

PRINT  1  LINE  THUS 

PR(BUSY  LINES)   PR( INPUT  LINES) 

SKIP  2  OUTPUT  LINES 

FOR  I  =  1  TO  12,  DO 

PRINT  1  LINE  WITH  1-1,  STATE. PROBLE( I) . 2  AND  QUEUE . STATS 

(I)/2  THUS 
**  *****  ***** 

•  • 

LOOP 

SKIP  5  OUTPUT  LINES 

PRINT  3  LINES  WITH  AVE. QUEUE. TIME .NO. ATTEMPTS  AND  NO. 

REQUESTS  THUS 
AVERAGEOPARS  QUEUE  TIME  =  **.**** 
NO.OPARS   ATTEMPTED   **** 
NO.OPARS   COMPLETED   **** 

FOR  I  =  1  TO  12,  DO 

LET  SUMMARY (I, RUN)  =  STATE . PROB( I)/2 

LOOP 

SCHEDULE  A  STATISTICS  IN  2880  MINUTES 

RESET  TOTALS  OF  LINES .BUSY, N. QUEUE  AND  QUEUE. TIME 

IF  RUN  EQUALS  10 

GO  TO  FINAL. STATS 

OTHERWISE 

ADD  1  TO  RUN 

LET  NO. ATTEMPTS  =  0 

LET  NO. REQUESTS   =  0 

RETURN 

'FINAL. STATS' 

START  NEW  PAGE 

PRINT  5  LINES  THUS 
BUSY  SYSTEM  SUMMARY 
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STATE  2/HR  4/HR  6/HR  8/HR  10/HR  12/HR  14/HR  16/HR  18/HR  20/HR 

FOR  I  =  1  TO  12,  DC 

PRINT  1  LINE  WITH  1-1,  SUMMARY  (1,1),  SUMMARY  (1,2), 

SUMMARY  (1,3) .SUMMARY (I, 4) ,SUMMARY( I , 5) ,SUMMARY( I , 6) , 

SUMMARY  (1,7),  SUMMARY(I,8) ,  SUMMARY  (1,9),  SUMMARY 

(1,10)  THUS 
**  .****  # ****   ****  < ****  ^****  #****  #****  #****   ******** 

LOOP 
STOP 

END 
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